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Abstract 

A new method for generating boundary-fitted, curvilinear, two-dimensional 
grids by the use of the Poisson equations is presented. Grids of C-type 
and 0-type have been made about airfoils and other shapes, with circular, 
rectangular, cascade-type, and other outer boundary shapes. Both viscous 
and inviscid spacings have been used. In all cases two important types 
of grid control can be exercised at both inner and outer boundaries. 

First is arbitrary control of the distances between the boundaries and 
the adjacent lines of the same coordinate family, i.e., "stand-off" dis- 
tances. Second is arbitrary control of the angles with which lines of 
the opposite coordinate family intersect the boundaries. Thus, both grid 
cell size (or aspect ratio) and grid cell skewness are controlled at 
boundaries. Reasonable cell size and shape are ensured even in cases 
wherein extreme boundary shapes would tend to cause skewness or poorly 
controlled grid spacing. An inherent feature of the Poisson equations 
is that lines in the interior of the grid smoothly connect the boundary 
points (the grid mapping functions are second-order differentiable). 

A user -oriented, well documented, FORTRAN computer program, called GRAPE, 
has been written to employ this grid generation method. It is available 
from the Applied Computational Aerodynamics Branch at NASA- Ames Research 
Center. 


*Now with Stanford University. 
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DESIRED PROPERTIES OF A GRID GENERATOR 




• ARBITRARY BOUNDARY SHAPES 

• ARBITRARY POINT DISTRIBUTION ON BOUNDARIES 

• SMOOTH VARIATION (DIFFERENTIABILITY) IN INTERIOR 

• COMPUTATIONALLY FAST 

• EASY TO USE 


• CONTROL OF ANGLES 
AT BOUNDARIES 


• CONTROL OF SPACING 
NEAR BOUNDARIES 



The principal contribution of this work 
is that the angles and spacing at the 
boundaries are input . Thus, one need 
not try to implement pre-determined 
angles and spacing by trial -and -error 
tuning of other parameters. 


TOPOLOGY OF GRID MAPPINGS 




PHYSICAL SPACE 


O-TYPE GRIDS 



PHYSICAL SPACE COMPUTATIONAL SPACE 

C-TYPE GRIDS 


Topology and notation for 0-type and 
C-type grids are shown here. The 
independent variables in the physical 
space are x and y, while £ and n are the 
independent variables in the cartesian 
computational space. 
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POISSON EQUATIONS 


£xx + £yy ^ 

V xx + *?yy * Q(t l) 

OR, WITH DEPENDENT AND INDEPENDENT VARIABLES INTERCHANGED: 

<**££ - 20x il} + yx vrj = -J 2 (PX| + Qx^J 
ay^ - 20yfy + yy^ = -J 2 (Py^ + Qy^) 

WHERE: a=x 7? 2 + y J? 2 

P = ^rj + V^r] 

7 7 

7 = x £ + 

J = x £ v » ‘ VS 


Basic to the method is that the grid transformations 
£ = £(x,y), n = n(x,y) must satisfy the Poisson 
equations. The equations are solved with dependent 
and independent variables interchanged to facilitate 
numerical integration and the application of boundary 
conditions. The equations with variables thusly 
interchanged are sometimes referred to as the "trans- 
formed" Poisson equations. 


CHOICE OF INHOMOGENEOUS TERMS 


P(£, v) - p(£)e” a7? + r(£ ) e^^max” 7 ^ 
Q(£, r?) = q(£) e~ b7? + s(£) e^^max” 7 ^ 


Inhomogeneous, or right-hand-side, or P and Q terms in 
the Poisson equations determine the character of the 
grid. Different choices for P and Q produce different 
grids. In this method P and Q are chosen as shown here 
with a,b,c, and d positive Note that the inner (n = 0) 
boundary PU,n) reduces to p(s), and that at the outer 
(n = omax) boundary P(c,n) becomes r(t), and similarly 
for QU>n). The approach is to assume that the geo- 
metric input requirements (control of angles and spac- 
ing at boundaries) are satisfied along with the Poisson 
equations at the boundaries, then back-sol ve for pUK 
qU) , rU), and s(e). Then PU,n) and QU,n) can be 
calculated for every point in the field. 


ITERATIVE UPDATE OF INHOMOGENEOUS TERMS 


GIVEN TWO GEOMETRIC REQUIREMENTS: 

• CONTROL 0(f) 

• CONTROL As(£) 


TWO ADDITIONAL EQUATIONS: 

’Vr? = I V£ I I V r) I cos 6 
As = [(Ax ) 2 + (Ay ) 2 ] 1/2 


' r 


NEW VALUES FOR p(£), q(f) 


\ 



DESIRED VALUES FOR x , y 
AT BOUNDARIES 


x vv , v vr) AT BOUNDARIES FROM 


NEW x, y 


NEW VALUES FOR P(£, r>), Q(f, t?) 



For each boundary (inner and outer) the two geometric control 
requirements can be re-cast as the two additional equations 
shown on the upper right. These two equations can be solved for 

the two derivatives x and y . Derivatives x , y , x , and y 

n n s s 

can be found by differencing known, fixed boundary data. Deriva- 
tives x.. and y^ are found by differencing the x and y , just 
Cn J y n n 

found, with respect to f. Thus, to back-sol ve the transformed 
Poisson equations for p(^) and q(f), two derivatives remain to be 

found: x and y . 
no nn 

In the solution procedure, each iteration step is in two parts. 
First, the x and y from the previous iteration (or initial condi- 
tions) are differenced to find x and y at the boundaries. These 

nn on 

are combined with all the other derivatives discussed above, which 
are fixed for all iteration levels, to form the transformed Poisson 
equations at the boundary. These are back-solved for new values of 
p(f) and q(f). Terms r(f) and s(f) are similarly found. New 
values for P(?,n) and Q(f,n) can then be calculated. The second 
part of each solution step is to perform one iteration of some 
solution procedure, such as SLOR. The above is iterated to 
convergence, producing a grid that satisfies the given geometric 
requirements. Inhomogeneous terms which yield the desired grid 
control are thus found automatically as the solution proceeds. 
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COARSE-FINE SEQUENCING SOLUTION PROCEDURE 




OF POOR Q'JALIT V 


Numerical convergence is greatly accelerated by an 
additional feature: coarse-fine sequencing. The 

solution is first iterated to convergence on a coarse 
grid consisting of every third point in the e, direc- 
tion and every third point in the n direction. This 
convergence requires relatively little computer time 
since the amount of arithmetic being done per step 
is one-ninth that which would otherwise be done. 

The coarse solution is then interpolated to provide 
initial conditions for a fine solution using all of 
the points. Coarse-fine sequencing produces a speed- 
up over normal SLOR by a factor of up to 15. Grids 
have been generated, for simple cases, in as little 
as two-thirds of a second of CPU time per thousand 
grid points on a CDC 7600 computer, including "set-up" 
overhead . 
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GRIDS ABOUT HIGHLY CAMBERED ELLIPSE 



GENERATED BY 
LAPLACE 
EQUATION, 
SHOWING POOR 
CONTROL OF CELL 
SIZE AND 
SKEWNESS 



GENERATED BY 
POISSON 
EQUATION 
WITH CONTROL 
AT BOUNDARIES 



The effectiveness of the grid control is demonstrated in 
this comparison of two grids about a highly cambered 12:1 
elliptical airfoil. In the two top figures is seen a grid 
generated by the Laplace equations — like the Poisson equa- 
tions but with P = Q = 0. Uncontrolled cell size and skew- 
ness are clearly seen. The figures on the bottom were 
generated by the present method with the grid control at 
the boundaries. It was required that the lines intersect 
the airfoil at 90° and that the standoff distance be 0.005 
chord lengths at all points on the airfoil surface. The 
angle requirement was satisfied to a tolerance of + 0.1° 
and the distance to a tolerance of + 0.00001 chord lengths. 






NEARLY SQUARE GRID CELLS ALONG SURFACE IN LEADING 
EDGE REGION OF NACA 0012 AIRFOIL 



GRID SPACING NORMAL TO SURFACE (As) SPECIFIED AS 
EQUAL TO LOCAL VALUE OF ARC-LENGTH 
ALONG SURFACE 


An interesting capability of this method is seen in this 
close-up view of the leading edge region of a grid about 
an NACA 0012 airfoil. The angle requirement need not be 
90° and it need not be equal at all points on the airfoil 
surface. Likewise, the standoff spacing need not be con- 
stant. In this grid the angle requirement was chosen as 
90° everywhere on the airfoil, but the standoff spacing 
requirement (normal to the airfoil surface) was taken to 
be equal to the local value of the arc-length al ong the 
surface. This produces grid cells along the surface that 
are nearly square, despite greatly varying surface dis- 
tribution. 



C-TYPE GRID FOR MODELING WIND TUNNEL WALLS 



GRID SPACING AND ANGLES CONTROLLED AT 
OUTER BOUNDARY 


A C-type grid for modeling flow through 
a wind tunnel is seen here. Grid spac- 
ing and angles are controlled at the 
outer boundary. 



IS 





UPSTREAM REGION IN O-TYPE CASCADE GRID 



CONTROL OF SPACING AND 
ANGLES AT OUTER BOUNDARY 
ENSURES SMOOTH TRANSITION 
BETWEEN CASCADE ELEMENTS 


Control of spacing and angles at outer 
boundary is applied here to a cascade. 

It was specified that the lines of con- 
stant n which intersect the top and 
bottom parts of the outer boundary of 
each cascade element do so vertically. 
Thus, the application of periodic bound- 
ary conditions between cascade elements 
is facilitated. 

This same capability ensures smooth ver- 
tical transition across the branch-cut 
in the wake region of a C-type grid. 
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GRID USED IN STUDY OF BLAST WAVE ENCOUNTERING 
STATIONARY PICKUP TRUCK 
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FEATURES OF GRAPE (GRIDS ABOUT AIRFOILS USING 
POISSON S EQUATION), A USER-ORIENTED 
FORTRAN COMPUTER PROGRAM 


• Q , As ARE INPUT 

• CODING IS MODULAR, GENEROUSLY COMMENTED, SYNTACTICALLY 
CONSERVATIVE 

• BUILT-IN DEFAULT CASE, AND SIMPLE, WELL THOUGHT-OUT INPUT WHICH 
IS CHECKED BEFORE USE 

• VERSATILE: C-TYPE OR O-TYPE, FREE-STREAM OR WIND TUNNEL OR 
CASCADE, VISCOUS OR INVISCID 

• GRAPHICAL OUTPUT 

• FAST 

• WELL DOCUMENTED AND ACTIVELY SUPPORTED 



A user -oriented , well documented, FORTRAN computer 
program, called GRAPE, has been written to employ 
this grid generation method. It is available from 
the Applied Computational Aerodynamics Branch at 
NASA-Ames Research Center. 



